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We study the temporal fluctuations in time-dependent stock prices (both individual and compos- 
ite) as a stochastic phenomenon using general techniques and methods of nonequilibrium statistical 
mechanics. In particular, we analyze stock price fluctuations as a non-Markovian stochastic process 
using the flrst-passage statistical concepts of persistence and survival. We report the results of 
empirical measurements of the normalized g-order correlation functions fq{t), survival probability 
S{t), and persistence probability P{t) for several stock market dynamical sets. We analyze both 
minute-to-minute and higher frequency stock market recordings (i.e., with the sampling time St of 
the order of days). We find that the fluctuating stock price is multifractal and the choice of St has no 
effect on the qualitative multifractal behavior displayed by the 1/g-dependence of the generalized 
Hurst exponent Hq associated with the power-law evolution of the correlation function fg (t) ~ t^'' . 
The probability S{t) of the stock price remaining above the average up to time t is very sensitive 
to the total measurement time tm and the sampling time. The probability P{t) of the stock not 
returning to the initial value within an interval t has a universal power-law behavior, P{t) ~ t~^ , 
with a persistence exponent 9 close to 0.5 that agrees with the prediction 9 = 1 — H2- The em- 
pirical financial stocks also present an interesting feature found in turbulent fiuids, the extended 
self-similarity. 
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I. INTRODUCTION 

The financial stocks are complex, nonlinear, open sys- 
tems characterized by a large number of parameters. 
Among other features, theypresent a multifractal be- 
havior [llllllllllalillllllSini Whether or 
not the multifractality is intrinsic or apparent is still 
an open question [l^ . Several interesting multifractal 
models '3, have been developed over the last decade. 
For example, the multifractal random walk model, in- 
troduced by Barcy and collaborators |3 has been 
recently shown to explain, besides the multifractality, 
other features of financial time-series, such as the ab- 
sence of correlations between price variations and long- 
range volatility correlations. Other multifractal models 
have been proposed to account for the asymmetry 

of the volatility-return correlation function. An impor- 
tant step toward a better understanding of such intri- 
cate dynamical processes is to search for new methods 
that are able to provide information about their tempo- 
ral evolution. One way to explore the temporal evolution 
of a stochastic system such as a fluctuating stock price, 
denoted by x{t), is by measuring the persistence prob- 
ability, P{t). That is the probability of the stochastic 
variable x{t) not reaching its original value correspond- 
ing to the starting time to up to a later time to + 1- This 
concept, closely related to the first-passage probability, 
has been successfully implemented in surface growth phe- 
nomena 15, 16] and has been used to determine the uni- 
versality class and the nonlinear features of the underly- 
ing dynamical process through the exponent 6 associated 
with the power-law decay P[t) ~ t^^ of the persistence 
probability at large times. Alternatively, one is inter- 
ested in measuring the survival probability S(t) which 



is defined as the probability of the stock price remain- 
ing above a reference value up to time t. In contrast 
with the persistence probability, wc show that the sur- 
vival probability depends independently on both the total 
measurement time tm and the time between successive 
recordings, i.e., the sampling time 5t. The concepts of 
persistence and first-passage time have been recently in- 
troduced in econophysics 0, 0| . Despite the capacity 
of first-passage statistical tools to predict the degree of 
performance of a given stock (i.e., how long a stock re- 
mains above a certain value, what is the first time when it 
reaches a particular level) , it is quite intriguing that their 
use in understanding the evolution of financial markets 
is rather scarce. Our study emphasizes the potential of 
the survival and persistence probabilities in any time se- 
ries investigations and shows how these concepts can be 
connected to the traditional analysis [3| based on the 
time evolution of price-price correlation functions. We 
also point out the extended self-similar behavior mani- 
fested by the stock price correlation functions. Extended 
self-similarity was originally observed in fluid turbulence 
problems 20, 21j and subsequently in discrete stochastic 
surface growth models [22|. Therefore we combine fea- 
tures from very different fields, such as surface growth, 
econophysics of stock markets, and fluid dynamics in an 
effort towards understanding the temporal evolution of 
flnancial stocks. 

Our statistical study of stock market temporal fluctu- 
ations is motivated primarily by the availability of the 
huge amount of quantitative data on the stock market 
prices both for individual stocks as well as for aggre- 
gate stock indices and both for instantaneous minute- 
to-minute price fluctuations as well as large-scale fluctu- 
ations over several years. Such detailed and precise quan- 
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titative information about the time dependence of a far- 
from-equilibrium stochastic process is not easy to find 
in real physical systems. For example, temporal thermal 
fluctuations of steps on solid surfaces, which we have re- 
cently analyzed (2^ 24, 25] to study the persistence and 
survival properties of equilibrium step fluctuations, usu- 
ally provide reliable time dependent data only over a cou- 
ple of decades. Similarly, experimental studies [2^ of ki- 
netic surface roughening in interface growth, which have 
served as model problems for developing the concepts 
of dynamical scaling in nonequilibrium growth phenom- 
ena, usually have reliable data for a couple of decades in 
growth time. The availability of extensive quantitative 
information on stock prices, including minute-to-minute 
price fluctuations within a single day as well as prices 
spanning over years (or even decades for some stocks and 
index averages), therefore provides an almost unique op- 
portunity to statistically analyze a time-dependent non- 
Markovian stochastic phenomenon using actual "experi- 
mental" data covering many decades in time, in contract 
to the corresponding experimental studies of real physical 
systems which span only a few decades in time. It should 
consequently be possible to obtain deflnitive information 
about the persistence and survival behavior of stock price 
fluctuations studied as a non-Markovian stochastic phe- 
nomenon. The availability of extensive stock price data 
should also allow for a detailed multifractal analysis of 
stock price fluctuations, which could serve, in principle, 
as a quantitative measure of "volatility" in the stock mar- 
ket since multifractality is directly related to intermit- 
tency and intermittency arises from the flnite probability 
for very large scale fluctuations. 

Our work presented in this article should therefore be 
thought of as a detailed empirical "experimental" anal- 
ysis of a non-Markovian dynamical stochastic process, 
namely the stock market, from the specific perspectives 
of first-passage statistics (i.e., persistence and survival 
probabilities and exponents) and intermittency (i.e., mul- 
tifractality and volatility). Our motivation for studying 
this specific problem is the existence of huge amount of 
experimental data (i.e., stock prices) of very wide dy- 
namical range. The fact that there is widespread current 
interest in the statistical physics community in the sub- 
ject of "econophysics" , which is nothing but the study 
of economics using the principles and methodologies of 
statistical mechanics, gives our current work some broad 
context, but our own interest is, however, complementary 
- we are using the vast amount of available dynamical 
data on stock price fiuctuations to carry out an "exper- 
imental" study of the important first-passage statistical 
concepts of persistence and survival in non-Markovian 
stochastic phenomena. The current work is, in some 
sense, a continuation of our earlier work on understand- 
ing various stochastic phenomena (i.e. thermal step fluc- 
tuations [2^ |3j us and nonequilibrium surface growth 
[1^ ) from the first-passage statistics perspective - here 
the stochastic process under consideration being an eco- 
nomic phenomenon (i.e. stock price fiuctuations) rather 



than physical phenomena as in the past. 



Our study of the multifractal character of price fluc- 
tuations is based on a multifractal versi on l27l| of the 
traditional detrended fluctuation analysis ^8|. We also 
use a standard dynamical scaling analysis, inspired from 
surface growth phenomena, to show the multifractal be- 
havior of the flnancial stocks, quantitatively extracted 
from the g-order price-price correlation functions. Our 
results go in parallel with earlier analyses of other groups 
QilMIE^IEIEQl which succeeded in showing multifrac- 
tality in several stock markets and commodities. On the 
other hand we investigate the persistence of price fluc- 
tuations described by the persistence exponent 9. The 
bridge between these two analyses is provided by the sec- 
ond order Hurst exponent H2 associated with the corre- 
lation function of the stock price, which has been shown 
[TEL to be simply related to the persistence exponent 
through H2 = 1 — 9. In this study we verify that this sim- 
ple relation is satisflcd by both low and high frequency 
fluctuating flnancial stocks. 



The data we use for our stochastic study comprise the 
daily Intel (INTO) stock value between January 1990 
and December 2002 and the composite NYSE index (i.e., 
NYA index) between January 1966 and December 2002. 
This corresponds to 3355 and 9312 data points, respec- 
tively. We also analyze sets of data recorded every minute 
for the Johnson and Johnson stock (JJ 2000) and ev- 
ery flve minutes for the Intel stock (INTO 2000) dur- 
ing the year of 2000 0. This corresponds to 98280 
and 19656 data points, respectively. The daily recorded 
stocks have an obvious exponential increase of their prices 
over several years. Therefore the exponential drift of the 
background can be subtracted from the stock price and 
we deflne a new stochastic variable x{t) = x{t) — Xb{t), 
where the stochastic variable Xb{t) associated with the 
evolution of the background depends on two parameters 
a and b, Xb{t) = aexp(M) [30j. We analyze the per- 
sistence probability of both x{t) and x{t). Changing the 
background subtraction within reasonable limits does not 
affect our statistical conclusions about persistence expo- 
nents and/or multifractality. 



The rest of the paper is organized as follows. In Sec.lTTI 
we define the various dynamical correlation functions and 
related statistical quantities, as well as the various expo- 
nents to be used throughout our statistical analyses - 
Sec. ^ is important in introducing the methodology of 
our analyses; in Sec. lIIII we present our extensive results 
and discussions; our concluding remarks are exposed in 
SecHvl 
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II. NONEQUILIBRIUM STATISTICAL 
MECHANICS TECHNIQUES 

A. Price— price correlation functions and extended 
self-similarity 

The generalized g-order price-price correlation func- 
tion is defined as 

G,{t)^{\xito+t)~x{to)n'/'^, (1) 

where x{t) is the stock price and the average is over all 
the initial times to. Gq(t) has a power-law behavior 

G,(^)^^^^ (2) 

which defines the exponent hierarchy Hq, also called 
the generalized Hurst exponent. The price evolution is 
multif racial if the exponent hierarchy Hq varies with 
q, otherwise is fractal (in the theory of surface dynami- 
cal scaling referred to as multi-affine and self-afhne, re- 
spectively). In particular, for g = 2, we recover the 
fractional Brownian motion case described by the well- 
known Hurst exponent, < H2 < 1. A simple way of 
assigning the presence of multifractality in a stochastic 
stock market is by looking at the multifractal spectra, 
Tq — qHq — 1. For fractals Tq depend linearly on q. A 
nonlinear behavior of Tq vs q is considered a manifesta- 
tion of multifractality. 

The temporal behavior in Eq. |2Jl of the generalized 
price-price correlation functions is analogous to the spa- 
tial behavior of the g-order structure functions of turbu- 
lent fluids or the generalized height difference correlation 
functions, C,(r,t) = (|/i(x-Fr,t) - /i(x, - r^'' 
(for small distances r), corresponding to surface growth 
models. The multiscaling behavior is revealed by the q- 
dependence of the scaling exponents ^q. The extended 
self-similarity (ESS) behavior in both turbulent fluids 
and surface roughening models refers to the enhance- 
ment of the scaling region once the log|Cm('')| is plotted 
against log| Cp,(r )|, where m and n are two different posi- 
tive integers |22| • In this study we show that for different 
financial stocks Gm{t) depends on G„(i) in a power-law 
fashion and, although the temporal scaling domain is not 
necessarily enhanced as in the case of turbulent fluids or 
roughening models, the ESS behavior is clearly identified 
and this interesting feature of the financial stocks can be 
further exploited to understand the associated multifrac- 
tal character. 



B. Persistence Exponent of Fractional Brownian 
Motion 

The fractional Brownian motion (fBm) is one of the 
simplest stochastic models that can be used to model 
financial stocks or any time series with long-range mem- 
ory. Before proceeding further, we summarize the def- 
inition and the known first-passage property of a fBm. 



We then describe the multifractal detrended fluctuation 
analysis (MF-DFA) method which is used to estimate the 
generalized Hurst exponent. 

A stochastic process ^{t) (with zero mean {^{t)) = 
0) is called a fBm if its two-time correlation function 
C{ti,t2) = ([?(ii) - ^(^2)]^) is (i) stationary, i.e., de- 
pends only on the time difference \t2 — ti\ and (ii) grows 
asymptotically as a power law jSlj 

C{ti,t2) ^\t2-ti\^", \t2-ti\»l. (3) 

The parameter < < 1 is called the Hurst exponent 
that characterizes the fBm and (• ■ • ) denotes the expec- 
tation value over all realizations of the process £,{t). In 
order to match the notation throughout the paper, let 
us call the Hurst exponent H2 instead of H. For the 
sake of completeness we also mention that, alternatively, 
a zero mean stochastic process ^{t) is called a fBm if its 
autocorrelation function has the following expression: 

adti,t2) = {atint2))-ti''''+ti"'^-\t2-t,r\ (4) 

The zero crossing properties of a fBm have been studied 
extensively in the past |3^|3l,|3- In particular, assum- 
ing that £,{t = 0) = 0, we are interested in the probability 
P{t) that a fBm does not cross zero up to time t (i.e., 
the persistence probability): 

P{t) = Prob {e(t') > 0, yt'<t}. (5) 

In terms of the stochastic stock price variable x{t), char- 
acterized by a particular value x{to) at the initial time 
to, the probability of remaining always above that value 
up to time to +t (i.e., positive persistence) reads 

P+{t) ^PToh{x{to+t') > xito), Vt' <t}, (6) 

and, similarly, the negative persistence probability reads 

P_(t) = Prob{a;(to +t') < x{to), ^ t' < t} . (7) 

These definitions can alternatively be reformulated in 
terms of the cumulative time series of the discretized 
log-returns. Let rj = hi[xj+Aj /xj] be the discrete set 
of log-returns, with j = 0,1,.... The sampling time 
Aj is the interval between successive measurements. We 
define the cumulative log-returns, Rj = J2i=o''^i- Since 
Rj = In [xj 1 2:0] , the definition of the positive persistence 
probability, for example, becomes: 

P+(A^) = Prob{i?„ > 0,V < n < iV}. (8) 

In several studies of linear surface growth models, 
characterized by identical positive and negative persis- 
tence probabilit ies, it has been shown that P(t) decays 
as a power-law [lali^ at large times, Pit) ~ t~^ , with 
the steady state persistence exponent d obeying the re- 
lation 

9=1-H2. (9) 
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We note that this relation holds for any zero mean pro- 
cess (not necessarily Gaussian js^ls^) that satisfies the 
requirements (i) and (ii) above. Both analytic arguments 
as well as numerical simulations supporting the relation 
@ have been presented previously in the literature in 
the context of fluctuating interfaces. In this study we 
investigate the behavior of P{t) at large times and its 
dependence on the sampling time for both x{t) and x{t) 
stochastic variables. 

The persistence probability can be generalized Is^ 
using the persistent large deviations probability, P(t,s), 
defined as the probability for the "average sign" ^av of 
the stock price fluctuation to remain above a certain pre- 
assigned value "s" up to time t: 

Pit, s) = Prob { 5av(t') > s, W < t }, (10) 

where 

S'av(i) = t"^ / dt' sign [x{to + t') - a;(to)]. (11) 
Jo 

Since Sav{t) G [—1,1], the probability P{t,s) is defined 
for — 1 < s < 1. For s = 1 we recover our earlier simple 
definition of persistence, while for s = — 1 the probability 
P{t, s = —1) is trivially equal to unity for all t. However, 
for the remaining values of the average sign parameter s, 
— 1 < s < 1, the generalization of the persistence prob- 
ability provides new information through the family of 
persistent large deviations exponents, 9i, associated with 
the power-law behavior, P{t,s) ~ t^^^^'^\ at large time 
scales. 



C. Survival Probability 

Perhaps of more practical interest in evaluating the 
temporal trend of a financial stock is the probability of 
the stock's price remaining above (below) a certain refer- 
ence value up to a later time + 1, given that its initial 
value at time was above (below) that reference level, 
i.e., the positive (negative) survival probability S±{t). 
Let us denote by S{t) the average between the positive 
and negative survival probabilities. This statistical quan- 
tity offers a better picture of the likelihood of a given 
stock having a positive evolving trend, for example, with 
respect to a preassigned reference value. The definition 
of this probability reads: 

S+{t) = Prob {x{to + t') > X, V i' < i} , (12) 

where x is the reference price. For simplicity we consider 
x to be the average price over the measurement time, 
t„n but in general it can take any value between the 
minimum and the maximum values of x{t) for all the 
discrete times t up to the final measurement time. We 
will show that S{t) depends independently on both 
and 5t and the scaling with t/5t appears only when 5t/tm 
is a constant. The same type of behavior has been found 



recently in experimental thermal fluctuations of surface 
steps on Ag(lll), screw dislocations on the facets of Pb 
crystallites and Al-terminated Si(lll) surfaces OE^- 

D. Multifractal Detrended Fluctuations Analysis 
(MF-DFA) 

MF-DFA is a reliable method for analyzing correlated 
time series. It is known to provide the accurate values of 
the generalized Hurst exponents even for time series with 
small length, while other similar methods, such as the 
Hurst Rescaled Range Analysis overestimate those 
values in the case of small size series [s^ . 

Let Xi, i — 1, . . . ,Tj be the stochastic price variable 
recorded at discrete times i. The final transaction time 
is denoted by Tf. We denote by rj the log-return price 
variable, Vj = hi{xj+i/xj), j = 1,...,T, where T = 
Tf — 1. We estimate the cumulative time series of the 
log-return price variables, 

i 

^« = E(^^-^)' * = l>--->?^: (13) 

where r = 1/2^ the average value of the log- 

returns. The time series X(i) is divided into Nr disjoint 
segments /„ (n = 1, . . . ,iV^) of equal size r. Obviously, 
Nt = [T/t]. For each segment we calculate the local 
trend using a linear least -squared fit IV ('t-, t) = an + bn t, 
where t € In and n = 1, . . . , Nt- The local time series of 
the cumulative log-returns is simply Xrin^t) = X[(n — 
1)t + t]. Therefore, the variance is given by 

I ^ 

F'(n,r) = -Y,iXr{n,t)-Yr{n,tyf. (14) 
t=i 

F{n,T) is called the fluctuation function. In order 
to avoid disregarding some data points X{i) when the 
length T of the time series is not a multiple of the time 
lag T, one has to repeat these steps starting from the 
opposite end of the interval. In that case, Xri^n^t) in 
Eq. (dll becomes equal to X[T - {n - Nr)T + t], for 
n — Nr + 1, . . . , 2Nr. By averaging over all the segments 
/„ we finally obtain the correlation function of order q, 

By construction, since we use a linear fit for simplic- 
ity, Fg(r) is defined for r > 3. The scaling form of the 
correlation function Fq{T) ~ r^' provides the family of 
generalized Hurst exponents, Hq. For reasons that will 
become clearer very shortly we also introduce the dimen- 
sionless fluctuation hmction, /(n, r), defined by 

(7 
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where <J = \J l/T'X]t=i('"t ^ '''Y the standard deviation 
of the log-returns during the interval T . Therefore, the 
dimensionless gth order correlation function becomes 

Obviously, fq obeys the same scaling relation as Fq, 

h{T)^CHj"\ (18) 

where Ch, is a constant independent of the time lag r. 
However, for q = 2 (which corresponds to the usual DFA 
procedure), the expression of this constant is known ex- 
actly 0, 

/ 2 1 2 y/^ 

where H2 is the Hurst exponent of the fBm. The time 
evolution of fq^r), along with the analytical result for the 
coefficient can be used to understand the dynamics 
and memory of financial stocks. 



III. RESULTS AND DISCUSSIONS 

We first discuss the results concerning the price-price 
correlation functions calculated using Eq. In Fig.Q] 
we present the results of Gq{t) for the high frequency 
stocks and in Fig. [3 for the low frequency stocks. It is 
obvious that since the log-log plots of Gq{t) vs. t do 
not exhibit linear behavior over the entire time range the 
associated Hurst exponent Hq varies with time and the 
scaling of the correlation functions suffers many transient 
regimes. A good power-law dependence appears for q — 
2. Although for other values of q the deviations from 
power-law become visible it is clear that Hq decreases 
with q. For illustration purpose we have fitted certain 
portions of these log-log plots to obtain a qualitative 
view of the dependence of Hq on as shown in Fig.|31 
In this figure we have used a large range of values for 
q (i.e., q = 1/10, 1/9, ...,1,2,. ..,10). We notice that Hq 
depends linearly on 1/g for both small and large values 
of \/q. The stock with the smallest sampling time of 
6t = 1 minute (JJ 2000) displays an increase of Hq at 
large values of 1/q, while for the rest of the stocks Hq 
has the tendency to saturate as 1/q increases. 

When using the MF-DFA method to calculate the cor- 
relation functions of order q we note that the power-law 
behavior of /^(t) vs. t extends over longer time periods 
making the extraction of the exponent Hq more reliable. 
The results are shown in Figs. 0] and |S1 This can be easily 
seen in the case of the daily NYA index and INTC stocks. 
The power-law is seen for more than three decades in 
Fig. while limited power-laws spanning two decades 
only are seen in Fig. |21 
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FIG. 1: (Color online) Log-log plot of the generalized price- 
price correlation function Gq (t) vs. t corresponding to minute- 
to-minute INTC stock (top) and JJ stock (bottom) for q = 
1/8, 1/4, 1/2, 1, 2, 4, 8 from bottom to top in each panel. 



As we have already mentioned, the value of Hq is sensi- 
tive to the fitted region of the log-log plot of /^(r) vs. r. 
This issue requires special treatment. Our strategy was 
to take advantage of the analytic result of Eq. (|19|) in 
order to identify the time window over which the sec- 
ond order Hurst exponent can be extracted correctly. 
This can be achieved by adjusting H2 such that the 
agreement between the empirical results for /2(t) using 
Eq. l(T7|) and the theoretical curve predicted by Eq. l(TH|) 
with q = 2 and Eq. H19|) becomes very good. This 
procedure is shown in Fig. |S| Once the time window 
which gives the best agreement for q ^ 2 is identi- 
fied, we use it to extract Hq for a large set of q val- 
ues {q = ±8, ±4, ±2, ±1, ±1/2, ±1/4, ±1/8). The re- 
sults for the g-dependence of Hq based on the MF-DFA 
method have been used to calculate the multifractal spec- 
tra, Tq — qHq — 1, shown in Fig. [7| We observe that 
all spectra deviate from a linear g-dependence, which is 
an obvious manifestation of the multifractality in these 
stocks. In the inset of Fig. [71 we also show the 1/g- 
dependence of Hq. At positive q, the qualitative trend 
of the results is the same as in Fig. 13 It is interest- 
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FIG. 2: (Color online) Log-log plot of the generalized price- 
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FIG. 3: (Color online) The generalized Hurst exponent Hq 
vs. 1/q for the four stocks discussed in the paper. Hq behaves 
linearly with 1/q at both small and large values of 1/q. For 
the daily stocks a linear least-square fit has been applied to 
the first two decades, while for the high frequency stocks, the 
fitted regions were 25 < i < 500 for the INTC 2000 stock and 
25 < t < 400 for the JJ 2000 stock. 



FIG. 4: (Color online) Log-log plot of the normalized fluctu- 
ation function /q(r) vs. time lag r corresponding to INTC 
2000 stock (top) and J J 2000 stock (bottom). The MF- 
DFA method has been used. The curves shown correspond to 
q = —8, —4, —2, —1/2, 1/2, 2, 4, 8 from bottom to top in each 
panel. 



ing to point out that the empirical set with the smallest 
sampling time, JJ 2000, which in the case of the stan- 
dard price-price correlation function analysis has shown 
an increasing trend of Hq at large values of 1/g, does not 
present this trend anymore, Hq saturating quickly as 1/q 
increases. Negative values of q are accessible within the 
MF-DFA analysis. For large negative values of 1/q the 
generalized Hurst exponents saturate rather fast, as in 
the case of large positive values of 1/q. 

The fact that Hq (and in particular H2) changes with 
time is a clear indication of the multifractal character. 
In this context we mention that the so-called multifrac- 
tional Brownian motion could be alternatively used to 
model this feature [s^- 

The analogy of the stock market fluctuations and fluid 
turbulence has already been pointed out in the litera- 
ture -"-t is known that in turbulent fluids the energy 
dissipation rate shows violent fluctuations. Similarly, as 
shown in Fig. |S1 we flnd that the temporal evolution of 
the local stock price differences, \x{t+6t) — x(t)\, presents 
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FIG. 5: (Color online) Log-log plot of the normalized fluctu- 
ation function /^(r) vs. r corresponding to the daily INTC 
stock (top) and NYA index (bottom). The MF-DFA method 
has been used to calculate fq{r). The curves shown corre- 
spond to —8, —4, —2, —1/2, 1/2, 2, 4, 8 from bottom to top in 
each panel. 
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FIG. 6: (Color online) The normalized fluctuation function 
/2 (r) as a function of time for the daily INTC stock during the 
period 1990-2002. The straight line represents the theoretical 
curve based on Eqs. II18II and 11191 for H2 = 0.47. 
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FIG. 7: (Color online) The multifractal spectra Tq vs. q for 
the stocks discussed in the paper. The deviation of all curves 
from the linear dependence is a signature of the multifractal 
behavior. The inset shows the 1/g-dependence of Hq, keeping 
the same symbol-stock correspondence as in the main figure. 



strong fluctuations which represent the signature of inter- 
mittency. In addition, the self-extending similarity fea- 
tures of intermittent fluid turbulence have been shown 
to exist in spatial height correlation functions of the ki- 
netic surface roughening models. We show in Figs. |H| 
and ini that the extended self-similarity exhibited by the 
structure functions in fluid turbulence also shows up in 
the temporal behavior of the financial stocks correlation 
functions. This observation offers a connection between 
three distinct physical problems, apparently without any 
intuitive connections: fluid turbulence, surface roughen- 
ing and financial stocks. 

For exemplification purpose we plot in Figs. 1^ andllUI 
log|/„i(r)| vs. log|/„(T)| for m = —1/2 and 2 and n = 1. 
From the linear behavior of these plots it is obvious that 
fm(T) ~ [/«(''')]""": where the expectation value for the 
exponents is amn = Hm/H„. We find that our em- 
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FIG. 8: The local stock price difference \x{t + St) — x{t)\ vs. 
the dimensionless variable t/5t for the INTC 2000 recorded 
with sampling time St = 5 minutes. 
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FIG. 9: (Color online) The normalized fluctuation functions 
/2(r) and /-i/2(-r) as a function of /i(r) for the INTC 2000 
and JJ 2000 (inset) stocks. 
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FIG. 10: (Color online) The normalized fluctuation functions 
/2(t) and f-i/2(j) as a function of /i(r) for the daily NYA 
and INTC (inset) stocks. 



pirical results for amn are in very good agreement with 
the expected ratios between Hm and as summarized 
in Table It would be interesting to check the exis- 
tence of the extended self-similarity in other financial 
stocks. We mention that the analysis of correlation func- 
tions based on the extended self-similarity is known to 
provide reliable values of the ratios between several gen- 
eralized Hurst exponents [2^ . 

Next, we present in the results of the persistence prob- 
abilities and persistence exponents for our fluctuating 
stocks. In Fig. we show the results based on the 
minute-to-minute stocks and in Fig.^]the daily stocks. 
We observe that the best power-law appears for the av- 
erage persistence probability P{t) = l/2(P+(t) + P-{t)), 
while departures from the power-law behavior can be 



Stock 


Hi /Hi 


H i/Hi 

■> 


a2,i 




INTC 


0.966 


1.051 


0.956 


1.065 


NYA 


0.960 


1.059 


0.961 


1.058 


INTC 2000 


0.958 


1.062 


0.956 


1.065 


JJ 2000 


0.983 


1.033 


0.962 


1.040 



TABLE I: Extended self-similar behavior of flnancial stocks 
based on the power-law dependence of log|/m(T)| on 
log|/„(r)| (m = — I and 2 and n = 1). The exponents a_i ^ 
and Q2,i extracted from the power-laws shown in Figs. Eland 
1101 are in very good agreement with the expected ratios of 
H_i/H\ and H2/H1, respectively. 



Stock 


e 


H2 


INTC 


0.51 


0.47 


NYA 


0.50 


0.50 


INTC 2000 


0.50 


0.47 


JJ 2000 


0.52 


0.47 



TABLE II: The persistence exponent 6 associated with the 
power-law decay of the average probability P{i) {P{t)) for 
the minute-to-minute stocks (daily stocks) . H2 is the second 
order Hurst exponent extracted from the time evolution of 
/2(t) (see Eq. GHJ). 



seen for P±{t) corresponding to the INTC 2000 stock 
and more clearly for the daily stocks. For the low 
frequency stocks, in addition to measuring the posi- 
tive, negative, and the average persistence probabilities 
of the stochastic price variables, we have also consid- 
ered the set of these three probabilities corresponding 
to the empirical sets after the background elimination, 
i.e., x(t) = x{t) — acxp{bt). For the INTC stock we 
have that a = 1.1703 and b = 1.1889 x 10"^^ and for 
NYA index a = 28.473 and b = 0.3177 x 10"^. The per- 
sistence curves for the variable x{t) are very similar in 
the sense that no distinction between the positive, nega- 
tive and average probabilities can be made. This result 
agrees with previous studies of the persistence probabil- 
ity of the German stock index We have used P{t) 
and P{t) in order to extract the persistence exponents 
for the minute-to-minute and daily stocks, respectively. 
The results are summarized in Tab. We compared 
these values against 1 ~ H2, with H2 extracted from the 
fitted power-law of f2{T), m order to investigate the va- 
lidity of Eq. ©. We find a good agreement between 6 
and 1 — i?2- However, it is important to emphasize that 
since both 6 and H2 are very close to 1/2 the memory 
effects of the time series under investigation can only be 
revealed by higher order correlation functions. The sec- 
ond order correlation function by itself cannot explain 
the multifractality discussed in this study since it indi- 
cates that the returns are uncorrelated. We also add that 
P{t) is not sensitive to the large discrepancy between St 
and tm corresponding to the high frequency stocks and 
daily stocks, respectively. 



9 



10 



10" 



I I I I I 1 — I — I — I I I I I I 1 — I — I — I I I I I 1 — n 



INTC 2000 



lL 



10' 



10^ 
t (minutes) 



10^ 





t (days) 




, p,-3 I I I I I I I I I I 

I 2 3 

10 10 10 10 

t (days) 



t (minutes) 



FIG. 11: (Color online) Persistence probabilities P{t) vs. t for 
the minute-to-minute stoclts. The dashed line corresponds to 
P+{t), the dotted line corresponds to P-{t), and the solid line 
represents P{t). 



The generalization of the persistence probability is 
shown in Fig. ^| We only present the results for the 
probability of persistent large deviations corresponding 
to the INTC 2000 stock, but we have checked the appli- 
cability of this concept to other empirical stocks as well. 
From the linear behavior of logP(f, s) vs. logi we con- 
clude that the ^-dependence of P{t, s) is indeed a power- 
law. We see that the local slope decreases as the average 
sign parameter s decreases. We have varied s from —1 to 
1 with an increment of 0.1 and the s-dependence of the 
resulting family of persistent large deviations exponents 
is shown in the inset of Fig. ^] We mention that each 
curve in Fig. 1131 corresponds to the average between the 
positive and negative persistent large deviations prob- 
abilities, i.e., Pit,s) = l/2{P+{t,s) + P^{t,s)). Both 
P+{t, s) and P-{t, s) show departures from the expected 
power-laws at large t, as we have seen in Figs. El and 
IT^ in the case of the positive and negative persistence 
probabilities. 

Finally we show the results of the survival probability 
in Fig. El We have looked at the temporal evolution 
of the JJ 2000 stock price recorded with two different 



FIG. 12: (Color online) Persistence probabilities P{t) vs. t for 
the daily stocks. Solid lines correspond to P+(t), P{t), and 
P-{t) (from top to bottom). The dashed lines correspond to 
P+ (t) and P- (t) . The circle represents the average probability 
P(t) for the variable x{t). 



10 LI I 1 1 1 1 — I — I I I I 1 1 1 1 — I — I I I I 1 1 1 1 — I — r-rr 

INTC 2000 ■ 




10 

t (minutes) 

FIG. 13: (Color online) Log-log plot of the persistent large 
deviations probability P{t, s) vs. t for the INTC 2000 stock. 
The average sign parameter, s, is varied from 1 to -0.8 (bottom 
to top) with an increment of 0.2. The inset contains the 
persistent large deviations exponent 9i vs. s. 
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sampling times, of 1 and 5 minutes, respectively. We 
have chosen different values of the measurement time, 
tm, which influences directly the value of the average 
price X — l/^m/Q™ dt' x{t'). As shown in upper panel of 
Fig. ll4l the survival probability measured over a longer tm 
has a slower decrease with time than the one correspond- 
ing to a smaller tm- We find that the empirical measure- 
ments of S{t) show good scaling with t/St at short times 
(t < 300 minutes), when the ration between the sam- 
pling time and measurement time is kept constant. This 
agrees with similar measurements done on experimental 
step fluctuations [23 | . Since a very simple interpretation 
of the stock market fluctuations is based on the random 
walk model, we have numerically simulated the survival 
behavior of a random walker which allows us to under- 
stand qualitatively all the features of S{t) found exper- 
imentally. Measurements of S{t) for the random walk 
model were carried out for systems of size L — 100. The 
measured average of the random walker variable at each 
site over the measurement time was used as the reference 
level in the calculation of the survival probability and the 
results were averaged over 300 independent runs. In bot- 
tom panel of Fig. we show that a perfect collapse of 
S{t) vs. t/5t appears when the ratio St/tm is constant. In 
addition, S{t) corresponding to the data recorded with 
the same St shows a slower decrease when the measure- 
ment time is larger, as in the empirical case. Therefore 
both sampling time and total measurement time have to 
be taken into consideration in order to interpret correctly 
the survival features of the financial systems. We want 
to point out that S(t) does not show an exponential be- 
havior over the investigated time range (see the inset of 
the bottom panel of Fig. I14|l . and possibly much larger 
t„i is needed to observe such a behavior at large time 
scales, as it happens in the case of equilibrium surface 
step fluctuations [25l |. 
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FIG. 14: (Color online) The average survival probability S{t) 
vs. t/5t for the JJ 2000 stock (top) with different sampling 
times and different measurement times, as shown in the leg- 
end. The bottom plot contains the numerical simulation of 
S{t) corresponding to a random walk with L = 100. The 
scaling of S(t) with t/5t appears when St/tm is constant. The 
inset shows the same curves on a logarithmic-linear scale. 



IV. CONCLUSIONS 

We conclude with some speculative thoughts on the 
possible development of "understanding" in the sense 
of physics with respect to stock price fluctuations. In 
physics, e.g. step fluctuations [ll El HHM El or 
kinetic surface roughening |l9l |. one typically looks for 
minimal (in the renormalization group sense) dynamical 
(in general, nonlinear) stochastic continuum partial dif- 
ferential equations underlying the stochastic phenomena, 
e.g., the Edwards-Wilkinson equation the MuUins- 
Herring equation j42] , the Kardar-Parisi-Zhan g equ ation 
[i^ . the Villain-Lai-Das Sarma equation [il. l45l| . etc., 
hoping to derive the long-time coarse-grained asymp- 
totic power law behavior of the system arising from 
some simple minimal underlying dynamics (which is of- 
ten based on symmetry and universality considerations). 
It is unclear whether such an approach based on con- 
tinuum nonlinear stochastic equations is at all meaning- 
ful for the understanding of the dynamical evolution of 



complex economic phenomena such as stock price fluc- 
tuations. It may be possible to empirically construct 
dynamical equations which are sufficient to reproduce 
the exponents and the related statistical stochastic be- 
havior described in this paper, but the necessary condi- 
tions for obtaining such equations are simply unknown 
(in fact, we do not know if such equations exist, except 
in some "trivial" data fitting sense). It is therefore quite 
intriguing that stock price fluctuation data are amenable 
to stochastic analyses based on first-passage statistics 
and multifractality (as carried out in this article) with 
results not that dissimilar from physical processes such 
as step fluctuations or nonequilibrium growth. 

Our work demonstrates that the persistence and the 
multifractal behaviors of stock prices (both individ- 
ual and composite) are subtle (including extended self- 
similar properties not identified in the literature before). 
It will be of interest to investigate if the empirical behav- 
ior we report in this paper can be derived from the various 
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multifractal stochastic models [Tol Ull ll^ IT^ for stock 
price fluctuations proposed in the literature. Such inves- 
tigations, clearly beyond the scope of the current work, 
would, however, be quite difficult since both persistence 
and extended self-similarity are notoriously difficult con- 
cepts to derive theoretically, even when the underlying 
non-Markovian dynamics is known for a process. At this 
early stage of our understanding of econophysics, the fact 
that the stock price fluctuations seem to follow the per- 
sistence and the multifractal properties of well-studied 
surface fluctuation phenomena is by itself intriguing and 
interesting. 

To summarize, in this study we have analyzed the mul- 
tifractality, extended self-similarity, and flrst-passage 
properties of several flnancial stocks. While the second 
order Hurst exponent and the persistence exponent char- 
acterizing the power-law decay of the average persistence 
probability are not able to explain the long-term corre- 
lations in the investigated price time series, higher or- 



der correlation functions reveal much richer information 
about the complicated dynamics of such systems. We 
have shown that the persistence exponent is in agree- 
ment with 1 — H2 and does not depend on the sampling 
time and measurement time. However, the survival prob- 
ability has a nontrivial dependence on both 5t and im, 
presenting scaling with t/5t only when the ration 5t/tm 
is a constant. The numerical simulations of persistence, 
survival, and extended self-similarity features using dis- 
crete models of financial stocks remain on interesting 
open problem. 

The authors wish to thank A. Christian Silva and Vic- 
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are indebted to Professor Michael E. Fisher for his proof 
reading of our manuscript and his positive comments. 
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